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A new formalism will be presented in order to study real time evolution of quantum systems at finite tem- 
perature. Probability distributions for time-correlated observables will be studied non-perturbatively and fully 
quantized. This works for various systems, including ones with tunneling phenomena. We have obtained good 
results with some computational methods which can be used on models with several degrees of freedom. Thus it 
looks feasible to study vacuum tunneling in real time for relevant field theories at finite T. 



1. Introduction 

We have been working on a non-pcrturbative 
approach to combine real time and finite Tem- 
perature for the full quantized theory. This is 
crucial for the study of systems with tunneling at 
intermediate temperatures. 

We choose to work directly with projectors onto 
initial and final states, given that the system at fi- 
nite T is described by the density matrix p. This 
method gives a direct probability interpretation 
for going from states where some degrees of free- 
dom have some values at time 0, to some other 
values at time t (in the presence of the thermal 
bath produced by all the degrees of freedom). 
This procedure is physically appealing as it multi- 
plies probabilities to be in quantum states (given 
by p) by the probabilities to go from these states 
(projected onto some states for the subsystem) to 
some states at later time. This contrasts with the 
usual < exp {-l3H) [Q(t) - Q(0)]2 > to estimate 
time correlations |l[] , where the Operator to mea- 
sure is positive by construction but one has the 
problem that it is hard to interpret as a transition 
probability. 

The probabilities to find the system in some 
configuration can be written in a simple way in 
terms of the path integral formulation in "com- 
plex" time. We have been able to solve several 
systems like the h.osc, the double well, (1 — cosx) 
potential or the interesting case of a ring of N 
coupled oscillators, with the aim of the field the- 
ory case with vacuum tunneling. 
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2. Finite time formalism 

Usually, at finite T = 1/(3, one just takes the 
trace over Euclidean time (3 to calculate expec- 
tation values. For time correlations we will need 
a path in the complex time plane of the shape: 
(0,0)^(t,-e)^(0,-2e)^ (0,-/3). 
We will combine finite T and real t non- 
perturbatively, by inserting suitable operators at 
the points (0,0), {t,-e) and (0, -2e). With a 
similar path Semenoff, Niemi, Weiss and others[P 
deduced Feynman perturbation rules. 

To fix ideas, let us ask for the probability that 
if the first particle (of a system of N interacting 
quantum particles) is located close to xi at t = 0, 
it will be close to j/i at time t, the system being 
at T = 1//3. We will see that this is given by: 

P{xi,yi)^ jdx-2..dx'2-dy2: Gtixi,X2, ■■;yi,y2, ■■) 

■ Gt'^{yi,y2, ■■;Xl,x'2, ..) Gf!{xi,X2, ■■;Xi,X2, ..)(1) 

In order to deduce this equation, we start from 
the knowledge that for t — the system is in a 
mixed state, described by the density matrix: 

p = exp(-/3H) = ^exp(-/3£;„)|n >< n\ (2) 



with the additional condition that the first degree 
of freedom be at xi, which is obtained by apply- 
ing the projector Px^ := \xi >< xi \ * ll2,..Af on 
both sides of p. So, for t ~ we describe the sys- 
tem by PxipPxi- This operator evolves as usual, 
U+{Px,pPx^)Ut , with Ut - exp(-ii/t) . 
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We measure again at a time t, with a position 
yi for the first particle. The probabiUty is then: 

Pfi^uyi) = ^^Tr[Py,U+{P.,pP.jUt] (3) 

where we have discarded one of the Py-^ due to 
the cychcity of the trace. Note that we have to 
take a smaU window around xi as otherwise we 
would measure the position with infinite precision 
implying total delocalization afterwards (we can 
see this effect numerically, in that the P{xi,yi) 
goes to a constant in yi for very small win- 
dows). This can be done easily by defining the 
projector on a window around the coordinate Xi: 
J^^^^ dxi\xi X Xi\ , which still has P^^ = P^^- 
The Normalization can be taken so that the inte- 
gration over yi in Eq.(3) is one, corresponding to 
the assumption that the particle is in the window 
[xi—a,xi+a], or with Z, in which case we include 
the a priori probability to be in that window. In 
both cases the integrated P is i-independent. 

Notice that classically Px^ can be commuted 
with p so that Eq. (3) would resemble a time cor- 
relation among projectors at finite (3, but in the 
quantum case this ordering gives a probability in- 
terpretation, while neglecting one of the P^^ one 
is mixing amplitudes for different states in the 
thermal mixture, which one should not do... 

We can recast Eq.(3) to show that it corre- 
sponds to the sum of probabilities to be in states 
|n > with position xi for one of the degrees of 
freedom by the probability to have gone to yi in 
time t : 

P{xi,yi)= j dyi^exp{-f3En) 

Ji=2,..,n „ 

• I / dxiGt(,yi,yi;xi,Xi)^n{xi,Xi)\^ (4) 

J i=2,..,n 

{Gt{x; y) :=< ?/|J7t|x > is the Green function and 
ipnix) a wave function) as we wanted to prove. 
Technically it is convenient to rewrite Eq.(3) as 

P(xi,yi) ^ / dxidx'^dyi < xi,Xi\p\xi,Xi > 

J i=2,..,n 

■ < xi,x'^\U^\yi,yi >< yi,yi\Ut\xi,Xi >, (5) 

a product of three Green functions as announced 
in Eq.(l). The point is that each of these Green 



functions can be calculated with the method of 
multiplying the transfer matrix for infinitesimal 
"dt" by itself n times to obtain times t = dt*(2"). 
This method works even for real t as we showed 
in our work with DuUweber and Hilf Q . Alterna- 
tively, one can implement a Fast Fourier Trans- 
form method Q to calculate the short time prop- 
agator and iterate for longer times. Both of these 
methods get improved nicely when introducing 
a small imaginary time ie in the time propaga- 
tion, as then the propagators are well behaved at 
large distances (suppressing the fast oscillations). 
With this ie we have even been able to use nor- 
mal Monte Carlo methods in real time, as now 
there is some weighting in the sampling, obtain- 
ing satisfactory results. The generalization to a 
field theory seems immediate, but technically it 
will be hard to implement even on coarse lattices 
due to the number of degrees of freedom. 

3. Solved models 

We have tested these methods for the systems: 

a) for a free particle, the probability to find the 
particle at some distance \yi—xi \ away has a near 
to gaussian shape, getting broader with time. As 
we higher the T, the dispersion grows faster. 

b) for the harmonic oscillator, as shown in fig.l, 
from the original window —.75 ±.5 for < = 0, it os- 
cillates with semiclassical frequency maintaining 
its coherence (as known for h.osc). With higher 
T the packet is broader but stays coherent. 

c) the double-well potential shown in fig. 2, 
where the original packet placed at one side dis- 
perses through tunneling to the other minima. At 
higher T the probability to be in the other minima 
increases, as expected. We are also studying the 
(1 — cos(x))" potentials with multiple minima to 
see the dispersion with t. This case is relevant, as 
similar tunneling via instantons in the field theory 
model is thought to produce the anomaly respon- 
sible for baryon asymmetry in the universe. 

d) The most interesting case that we have been 
able to solve, is the N particle chain of coupled 
oscillators in a circle, with transversal oscillation. 
This is the first genuine finite temperature system 
that we treated in the sense that the Temperature 
can be dynamically generated by the system in a 
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self consistent way (in contrast with an external 
given heat bath). 

We have been able to solve this model " almost" 
analytically by going to normal modes, Xi ^ X^- 
In these new coordinates the Hamiltonian is sep- 
arable and the N-particle Green function is just 

TV 

< x,,.,XN;t\x[, .,x'j^;0 >=n< ^.;i|^^;0 >. 

The transformation being linear and with Jaco- 

bian 1, the Green's function in the original Xi is 
simply equal to the new one, with Xg = Xs{xi). 

As the integrations in Eq.(l) have to be done 
numerically we have been able to solve for up to 
A'^ = 4 particles. The probability distribution is 
again some Gaussian-like curve that gets broader 
with time (the zero mode corresponding to a col- 
lective movement of the circle of strings is allowed 
and so the center of mass of the system disperses 
as usual in QM). As we increase the tempera- 
ture the coordinate xi disperses much faster due 
to the momentum transfer from the other oscil- 
lators. It would be interesting in the future to 
study to which extent we have reached the ther- 
modynamic limit in the sense that if we have a 
very large circle of N oscillators the probability 
to find the first in some position after some time, 
should converge to some limit which should not 
change when we go to A'" + 1 oscillators. 
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Figure 1. Probability to be at y for an h.osc. for 
times t « 0, 7r/2 and tt, at low and high T. 
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Figure 2. Probability to be at y for a double well 
for f w 0, .3 and 1. of tunneling time, for two T. 

4. Conclusions 

We have been able to combine non perturba- 

tively real time and finite Temperature in first 
quantization. It can be seen as the product of 
three Green functions in the complex time plane. 
The "matrix iteration" and the FFT methods 
work well in calculating the real time path in- 
tegral, improving convergence with a small ie. 

There are several interesting problems that 
seem shortly doable, like the ammonia molecule 
or several spin models (easy due to discrete field 
values). Our ultimate aim is to find a way to 
treat the field theory case of tunneling with in- 
stantons in real time and finite T, which is in 
principle doable with these methods but will be 
hard to implement due to the number of degrees 
of freedom. 

We thank L. Policy for useful discussions and 
the RRZN in Hannover for supercomputer time. 
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